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Abstract 

We investigate, in terms of production from pulsars and their nebulae, the cosmic ray 
positron and electron fluxes above ~ 10 GeV, observed by the AMS-02 experiment up 
to 1 TeV. We concentrate on the Vela-X case. Starting from the gamma-ray photon 
spectrum of the source, generated via synchrotron and inverse Compton processes, we 
estimated the electron and positron injection spectra. Several features are fixed from 
observations of Vela-X and unknown parameters are borrowed from the Crab nebula. 
The particle spectra produced in the pulsar wind nebula are then propagated up to the 
Solar System, using a diffusion model. Differently from previous works, the omnidi¬ 
rectional intensity excess for electrons and positrons is obtained as a difference between 
the AMS-02 data and the corresponding local interstellar spectrum. An equal amount 
of electron and positron excess is observed and we interpreted this excess (above ~100 
GeV in the AMS-02 data) as a supply coming from Vela-X. The particle contribution is 
consistent with models predicting the gamma-ray emission at the source. The input of 
a few more young pulsars is also allowed, while below ~100 GeV more aged pulsars 
could be the main contributors. 
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1. Introduction 


The AMS-02 experiment extended the observed cosmic ray (CR) electron, positron 
and electron plus positron spectra from 0.5 GeV up to 700 GeV, 500 GeV and 1 TeV 
respectively [see 1,2]. CR particles generated and accelerated at known sources are 
considered as primaries. For instance, the main component of electron spectrum is 
that produced by supernova remnants (SNR). CRs are also produced directly inside the 
interstellar medium (ISM). In fact, positrons were supposed to be mainly originated 
from the decay of muons produced by CR interactions with the ISM [e.g., 3]. These 
particles are commonly referred to as secondaries. Primary plus secondary CR spectra 
outside the region interested by the solar activity (i.e., the heliosphere) are known as 
local interstellar spectra (LIS, see e.g.. Sect. 2). In this work we focus our attention on 
electron and positron spectra. Moreover, we will refer to electrons produced in SNR 
and in the ISM as the “classical” electron LIS and to positrons produced in the ISM as 
the “classical” positron LIS (e.g., “classical” LIS, hereafter cLIS). 

At low energy, less than ~ 10 GeV, due to solar modulation, the observed CRs 
spectra deviate from LIS’s (see Fig. 1 and, for instance, 4, 5). At higher energy, it is 
commonly acknowledged that, inside the heliosphere, particle propagation is little af¬ 
fected by solar modulation, thus the omnidirectional distribution is the one determined 
by the LIS. Nevertheless observed spectra of electrons and positrons [see 1, 2] exceed 
the cLIS computed with GALPROP (see Sect. 2) at high energies (e.g., see Fig. 1). We 
evaluated the excess for electrons and positrons subtracting the cLIS from the AMS-02 
fluxes and we found an equal amount for the electron excess and the positron one (see 
Fig. 2). Therefore, we take the electron plus positron flux observed by AMS-02 as 
reference data, extending the comparison with models up to 1 TeV due to the experi¬ 
mental accuracy. 

In this paper, we investigate possible astrophysical sources of positrons and elec¬ 
trons i.e., pulsars and their nebulae (see Sect. 3 and e.g., 6, 7, 8), which may account 
for the flux excess, without the need to look for more exotic explanations, e.g., in 
the framework of dark matter scenarios [see e.g., 8, 9, 10]. Previous works already 
explored this scenario (see e.g., 11, 12, 13, 14, 15, 16). Differently from the usual ap- 
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proaches, we evaluated the contribution from astrophysical sources reproducing photon 
spectra in agreement with the observations. 

In particular we consider Vela-X as a source of pair production and acceleration. 
The positron and electron injection spectra are obtained using the diffusion process in¬ 
side the pulsar wind nebula (PWN) of Vela-X (see Sect. 3.1 and e.g., 17). Shape and 
luminosity of the flux are fitted in order to be consistent with the gamma-rays spectrum 
observed from the source. Using a diffusion model described in [18], we evaluated the 
particle spectra at the Earth position (see Sect. 3.2). A comparison with experimental 
data is, finally, discussed in Sect. 4, while we report our considerations and conclusions 
in Sect. 5. 

Preliminary results, on materials obtained in this paper, were presented as confer¬ 
ence contributions in [19,20]. 

2. Electron and positron Spectra 

2.7. The “classical" Local Interstellar Spectra 

The propagation of cosmic rays in the Galaxy can be described by the diffusion 
equation (e.g., 21, Chap. 3 and 22): 

dfi d 

+ V • [a- Vn,] + — [hi ni\ - pini + Pi, (1) 

where the time evolution of the energy density n, = dNijdE of cosmic ray species i with 
energy E depends on the source term Qi, diffusion coefficient D, usually described by 
a power law in the energy D{E) - Do(EjE^f, [see 23, 24], the change of the particle 
energy per unit time 7>,, catastrophic processes p, and nuclei collisions P,. Equation 
(1) accounts for i) the propagation of primary components like, e.g., electrons, protons 
and carbon nuclei mainly accelerated in SNRs [21, Chap. 4] and ii) the production of 
secondary spectra like, e.g., positrons and Boron nuclei produced from interaction of 
primary CRs with the ISM. 

The most recent data provided by PAMELA and AMS-02 were discussed by many 
authors. Eor instance, [25] introduced two slopes in the diffusion coefficient in differ¬ 
ent energy range to explain the proton and helium spectrum; [26] interpreted the CR 
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spectra using an inhomogeneity source distribution following the stars concentration 
in the galactic spiral arms; [27] modelled the rising of the PAMELA and AMS-02 
positron flux as due to purely secondary origin, without taking into account the en¬ 
ergy losses. The GALPROP model solves numerically Eq. (1) for all the relevant CR 
species in a cylindrically symmetric space [28], with a galactic radius Rcai. and height 
hcai. ■ Therefore, hereafter, we will use the most comprehensive propagation model; 
i.e., the GALPROP model by which we evaluate, at the same time and with the same 
propagation parameters, the local interstellar spectra of several kind of particles: pro¬ 
tons, electrons, ions, anti-particles and photons used for the electron energy loss. 

The GALPROP code' returns the predicted “classical” LIS for the specific particle 
at the Solar System. The solution of Eq. (1) depends on parameters like the boundary 
conditions of the galactic effective volume for CRs diffusion, the diffusion coefhcient 
and the injection spectra characterized by power laws with different spectral indices 
for nuclei, protons (jp) and primary electrons (yg). To determine these parameters, 
we compared the spectra obtained in this way with the experimental data above ~ 10 
GeV (energies high enough to neglect solar modulation effects), then we tuned the co¬ 
efficients minimizing the discrepancies. The calculated spectra were normalized at 50 
GeV with measured proton, electron and ion fluxes at Earth. Lor proton and electron 
spectra we used the AMS-02 data [1, 29], while for the ions ratios B/C, Be/B, Be/C, 
Li/B, Li/Be and Li/C we referred to the online cosmic ray database reported in [30]. 

The available data are best described using the parameters listed in Tab. 1. In 
Lig. 1 we reported the comparison between the cLIS’s and AMS-02 data for electron, 
positron and electron plus positron spectra. The cLIS’s, coming from GALPROP, were 
reported, in solid lines, for energy above 10 GeV where solar modulation is negligible. 

Since the choice of the values of Tab. 1 is not unique, we modified, one by 
one, the main GALPROP parameters responsible for the diffused spectra (the galactic 
height and the diffusion coefficient), as reported in Tab. 2. The ranges of this param¬ 
eters were determined keeping the produced LIS within the experimental data errors, 
with special regard to the Boron over Carbon ratio [see e.g., 31, 32, 33, 34, 35]. This 


'http://galprop.stanford.edu/webrun.php 
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Figure 1: Electron, positron (top) and electron plus positron (bottom) omnidirectional intiensities observed 
by AMS-02 [1,2] and the cLIS’s evaluated using GALPROP (obtained with the parameters reported in Tab. 
1) in solid lines; the allowed range of the parameters (reported in Tab. 2) is kept into account by the shadowed 
bands. 







Parameter 

Value 

Real. 

30 kpc 

hcal. 

+4 kpc 

Do 

5.8 • 10^* cm^ s“' 

6 

0.33 

Eq 

4 GeV 

Va 

30 km s^' 

Jp 

1.98 (£ < 9 GeV), 2.42 (£ > 9 GeV) 

Je 

1.7 (£ < 4 GeV), 2.68 (£ > 4 GeV) 


Table 1: Propagation parameters used in GALPROP code to determine the “classical” electron and positron 
US’s. 


overall uncertainty is included in the shadowed regions of Fig. 1. We also found that 
the cLIS, obtained with the parameters of Tab. 1 and 2, is compatible with the pjp 
ratio of PAMELA [36] and AMS-02 [35]. 


Parameters 

Range 

Galactic height (kpc) 

Diffusion Coefficient Constant (cm^s ') 

Diffusion Coefficient Index 

2 < hcai. < 6 

4 ■ 10^*^ < Do < lO^*^ 

0.3 < (5 < 0.4 


Table 2: Ranges of propagation parameters used in GALPROP code to determine the errors in the LIS 
evaluation. 

2.2. Electron and positron flux excess at high energy 

The omnidirectional intensity excess for electrons and positrons are shown in Fig.2. 

The difference between the observed AMS-02 spectra and GALPROP cLIS’s (solid 
lines of Fig. 1) were calculated for energy above ~ 10 GeV (where the solar modu¬ 
lation effects are negligible) and requiring at least a difference (above 10%) between 
the two fluxes. Under these constraints the electron and positron signal is reported for 
energy above 90 GeV and above 10 GeV, respectively. We report also the electron 
plus positron spectrum, above 50 GeV, divided by a factor two for a comparison with 
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Figure 2: Omnidirectional intensity excess for electrons, positrons and half of all electron, obtained as a 
difference between the AMS-02 flux and the con'esponding “classical” LIS. 

respect to the other data. The error bars of these data come from the experimental ob¬ 
servations. We can remark how these excess spectra of positrons and electrons can be 
fitted using similar power laws. The electron signal spectral index, resulting from the 
fit, is -(2.503 + 0.353), for positrons we have -(2.502 + 0.030), while for electron plus 
positron spectrum we have -(2.568 + 0.088). The points of Fig. 2 are dependent on the 
parameters used in GALPROP. The uncertainties due to the choice of the GALPROP 
parameters result as a scale factor of the omnidirectional intensities in Fig. 2. This 
uncertainty is mostly constrained by the positron spectrum and can be accounted as a 
scale factor of ~ 5% at 100 GeV and above, while at lower energy it is < 20%. For 
the electron plus positron spectrum, the uncertainty at 1 TeV is about 25%. Hereafter, 
we will compare our source models with the electron plus positron spectrum since the 
upper point reaches 1 TeV. 
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3. Pulsars as possible sources of the excess components 


Looking at Fig. 2, above 50 GeV, a region marginally affected by the solar modu¬ 
lation, it is possible to say that the electron and positron signals are compatible. Thus, 
the dominant physical process is expected to proceed via a pair production mechanism. 
Pulsars (PSRs) are among the most likely sources of electron-positron pairs. The high 
magnetic field and the fast rotation of these neutron stars lead to huge electric fields 
[37]. Electrons produce high energy photons via curvature radiation. The interactions 
between photons and high magnetic fields generate pairs [38] that can generate again 
curvature photons. At the end, electromagnetic showers are produced. 

PWN identify the region around the pulsar where a relativistic magnetized wind 
is populated with electrons and positrons [e.g., 39, 40]. PWN are widely believed 
to be responsible for the acceleration of cosmic rays up to energies of 10^^ eV [e.g., 
41, 42, 43]. The central pulsar converts its spin-down power into a relativistic wind 
injected near the magnetosphere. The electrons and positrons in the wind, interacting 
with shock fronts, are accelerated and get a power-law energy spectrum. They radiate 
at lower energies, from radio frequencies to X-rays, through the synchrotron process in 
the magnetic field of the nebula. The highest energy part of the spectrum comes from 
inverse Compton scattering of the radiation field; synchrotron radiation generated by 
the PWN, cosmic microwave background (CMB), infrared and star-light photons (e.g., 
44, 45, 46). The gamma ray spectra of PWN reach very high energy, for instance, tens 
of TeV. The observed high-energy emission from the Crab Nebula has been modelled in 
detail by several authors [see 44, 45, 47, 46, 48]. High-energy processes in other PWN 
such as that of the Vela-X nebula, the nebulae around PSR 1706-44, PSR 1509-58, 3C 
58, CTB 80 and other few nebulae have been also studied in [49, 50, 51, 48, 52]. It was 
also suggested that the production of gamma-rays in the interactions of hadrons with 
the matter of the supernova could contribute to the higher energy end of the observed 
spectrum, especially in the case of the youngest nebulae [see 53, 45, 54, 52]. 



3.1. Electron and positron spectra at the source 

Since pulsars are rotating magnetized neutron stars, basically PSR lose their en¬ 
ergy through the electromagnetic radiation of the magnetic dipole spinning around a 
tilted axis. Within a simplified model, it is possible to determine the main physical 
parameters of the pulsar (like magnetic field or energy loss) starting from the rotation 
frequency (v) of the pulsar, its first (v) and second (v) derivative. The spin down of the 
pulsar is assumed as: 

V = -kv", (2) 


where k depends on the magnetic moment and on the moment of inertia of the neutron 
star and n is the braking index that would be equal to 3 for a magnetic dipole. Assuming 
constant k, from Eq. (2) we get: 


,V2- 


( 3 ) 


Moreover, integrating Eq. (2), assuming constant k, the age of the pulsar (Tags) can be 
obtained: 


‘age 


1 - 1 - 

vo 


( 4 ) 


(n - l)v 

Here vq is the birth rotation frequency. The time evolution of spin-down luminosity is 
given by [e.g., 55] 

-U+i 

L(f) = To(l + -'l ( 5 ) 


\ TO/ 

where Lq is the initial spin-down luminosity and tq = Po/(n - l)Po is the initial spin- 
down time scale. 

To evaluate the energy spectrum of electrons and positrons at the source, we follow the 
approach described in [17]. This method is similar with others reported in [56, 55]. 
The time dependent model by [17] describes both particle and photon injection spectra 
from the nebula. It assumes an initial particle injection rate, produced by the pulsar, 
that follows a broken power law with indices ai and a 2 and energy break Et,: 


Q(Ee,t) 


QaiOiEJEhT' if E,< El, 
Qo(t)(EJEhr if Ee> Eh 


( 6 ) 
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where Eg is the particle kinetic energy and Qoit) can be derived by requiring the conti¬ 
nuity of the two power laws and that 


'' 


Q{Ee, t)EgdEg = 


(7) 


with rj the conversion factor of the spin-down power L{t) into particle luminosity. The 
diffusion equation for the differential electron density rigiEg, f) - dN(Ee, t)ldEe can be 
approximated as ( [17, 56, 55]): 


dne(Eg,t) ng(Eg,f) ng{Eg,t) 

- - -= Q{Ee, t) - 

dt 


(8) 


'^esc it) 

The particle spectrum obtained as a solution of the Eq. (8) over time from f = 0 to 
t — T (age of the PWN) is: 


dN(Eg, T) 
dEe 


f 


Q(Ee, t) exp 


T-t 

Teff 


dt. 


(9) 


where corresponding to the lifetime of an electron with respect to 

both synchrotron energy loss and escape timescale (i.e., the time to diffuse 1 PWN 
radius). In Ref. [17], the authors applied their model to the Crab nebula, for which 
an accurate spectral energy distribution is known, due to numerous observations from 
radio frequency up to gamma rays (see 17 and references therein). The low energy 
spectral index {a\ in eq.6) is generally fixed by fitting the synchrotron spectrum at low 
frequency. The high energy spectral index (q '2 in eq-6) is however related to the inverse 
Compton process. In the next subsection we will use this result as source spectrum of 
the electron and positron excess spectra at Earth. 


3.2. Electron and positron spectra at the Earth 

Eollowing the approach reported in [18], the time evolution of the energy density 
ng{x, £, f) of electrons or positrons from a single source distant x from the Solar System, 
with energy E and after a diffusion time f, is obtained disregarding the last two terms 
of Eq. (1) [see 57], i.e., 

dng(x, E, t) 


dt 


= Qe(E, f) + V ■ [oiEWngU E, f)] 


[KE)ne(x,E,t)]. 
oE 


( 10 ) 
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In Eq. (10), the term b{E) accounts for the rate of energy loss due to ionization, 
Bremsstrahlung, synchrotron and inverse Compton processes [e.g., 58, Chap. 4]. How¬ 
ever, above ~ 1 GeV, the only relevant mechanisms are synchrotron and inverse Comp¬ 
ton. Furthermore, above few GeV, using an average interstellar magnetic field of 3 /rG 
and the photon radiation fields reported in [57] (Table 2, model Ml), the fit of the total 
energy loss rate can be described by a power law as in; 

dE , 

_b(E)--boE\ ( 11 ) 

dt 

where bo ~ 1 ■ 10“'^ GeV^^s"^ (value in agreement with those reported in e.g., 59, 60). 
Due to the high rate of energy loss, a positron or an electron of 100 GeV dissipates most 
of its energy in about 10^ years and can diffusively travel up to a typical distance of 
about 2 kpc. Thus, sources responsible of the high energy positron and electron excess 
(observed by PAMELA, 61, and AMS-02, 62, 63) are located in a region relatively 
close to Earth (within a distance of ~ 2 kpc). 

Then, we fit the injection spectrum, i.e., the results of Eq. (9), as a function of the 
initial energy Eq, using a power law with spectral index a and an exponential energy 
cut-off Ecut'- 

e,(£o,f) = &.o(OVexp(-|^), (12) 

\ i^cut I 

the interstellar diffused spectra of electrons and positrons from Eq. (10) is [22, 18]: 



where t is the diffusion time (i.e., the time spent to reach the solar system), x is the 
distance between the source position and the Earth and Ad is the mean distance travelled 
by particles with initial energy Eq - EI{\-botE) down to energy E resulting from both 
energy loss and diffusion processes given by 


UE,Eo) 


I D(E')dE'\‘^^ 
[Je b(E') I 


(14) 
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Eq. (13) allows one to evaluate the electron and positron spectra at Earth coming 
from PWN for energies beyond ~ 10 GeV, i.e., energies high enough to neglect solar 
modulation effects. 

4. The case of Vela-X: Results 

The pulsar J0534+2200, located inside the Crab Nebula, the remnant of a supernova 
explosion occurring in A.D. 1054, is an extremely well studied object. This young 
source can give us information about the first step of the life of a generic PSR/PWN. 
Note that the Crab distance is about 2 kpc and, due to the age, we can not see yet 
particles coming from that source at Earth. The TeVCat catalogue^ contains less than 
40 PWN observed in the TeV energy range. Only five of them are closer than 2 kpc and 
were observed by Cherenkov telescope experiments like HESS [64, 65], Veritas [66, 
67] and Magic [68]. Vela-X belongs to this sample. These observations regard a small 
fraction of known pulsars and they are much less complete and accurate in comparison 
with the Crab Spectral Energy Distribution (SED). Indeed it is widely believed that 
PWN are no more observable after the early phase of expansion. [18] suggested that 
all the pulsars have an initial stage as PWN and the lifetime of these objects is about 
10^-10“* years. During this phase electrons and positrons are trapped inside the PWN, 
but later, after a time T from the SN explosion, they are free to propagate. Mature 
pulsars, like Geminga and Monogem, have no more gamma-ray emission from the 
nebula, but the electrons and positrons released are still coming to the Earth. Eor all the 
older pulsars we do not have information regarding the nebula photon spectrum, the 
braking index or the birth frequency. Eor what concerns the PSR age, we can roughly 
estimate the minimum characteristic age Tage,c - -v/2v as reported in [69]. 

4.1. Model 1: using observed parameters for Vela-X 

Vela-X was detected by HESS [65] in the very high energy gamma ray band, the 
spectrum can be fitted by a power law with the photon index Fy - 1.45 + 0.09+ 
Q.2sys in the energy range between 550 GeV and 65 TeV and an exponential cut-off 


^http://tevcat.uchicago.edu/ 
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at an energy of 13.8 + 23stat ± TeV. The X-ray part was detected using ROSAT 
combined with ASCA data [70, 71]. The spectrum observed in this region has a spectral 
index of Tx ~ 2. The observation of the timing property of the pulsar gives us the basic 
information to evaluate the braking index as in Eq. (3). Vela is one of the few pulsars 
[72] for which it is possible to measure the braking index: n - 1.4+0.2 [73]. Therefore, 
we built a model based on the observed parameters of Vela and we assumed a value of 
n not dependent on time. The birth period Pq can be evaluated by Eq. (4) knowing the 
frequency, its derivatives and the age of the object. The CRAB source is the only one 
for which all these parameters are known. Eor this reason, we assumed the same initial 
rotation period of the Crab pulsar (Pq ~ 20 ms, 74). In this way, Eq. (4) gives an age 
for Vela of about 26 kyrs (instead of the common characteristic age Tage.c = 11 kyrs). 
The results reported in Eig. 3 and 4 (Model 1) come from this analysis. We get a photon 
spectrum compatible with the HESS and ASCA data requiring a conversion efficiency 
of the spin down luminosity of about rj - 0.5% for both electrons and positrons. The 
Model 1, reported in the figures, is the diffused spectrum at the Earth for which we set 
r ~ 10 kyrs, in comparison with the initial spin-down time scale which is evaluated to 
be To ~ 29 kyrs. We also decided to vary the Po-value in the range: 10 < Po (ms) < 30. 
Consequently, Vela age changes from 33 to 19 kyrs. In all cases, we fit the photon 
spectrum in agreement with the ASCA and HESS data. Eor this purpose we change the 
efficiency inside the range 0.25 < rj (%) < I respectively and the released time from 
22 to 5 kyrs. The band of the model reflects the uncertainty on the Vela pulsar distance 
that is about 6% (Vela distance is 2%1"^^ pc) [69] and the variation due to the assumed 
initial rotation period 10 < Po (ms) < 30. 

4.2. Model 2: using Crab-like parameters for Vela-X 

Observation of the timing properties of the other five pulsars, younger than Vela, 
gives for the braking index n values in the range between 2 and 3 [72]. All of them are 
more similar to the value for the Crab nebula with respect to Vela-X. Therefore, we can 
alternatively assume that all the pulsars are similar at their birth and then we can take 
the properties (initial rotation period, breaking index) of Crab. We need to assume that 
there is a variation of the braking index from 2.5, at the birth, down to 1.4, at later time. 
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Parameter 

Model 1 

Model 2 

a\ 

1.9 

1.5 

ai 

2.8 

2.4 

Eh (MeV) 

1.5 • 10^ 

1.5 ■ 10^ 

n 

1.4 

2.5 

ri(%) 

0.25 - 1 

5 

Po (ms) 

10-30 

20 

Tage (kyr) 

33 - 19 

11 

To (kyr) 

22-35 

0.7 

T (kyr) 

5-22 

1 


Table 3: Parameters used in Model 1 and 2 for Vela-X nebula. 


It could be associated to some changes in the structure of the neutron star, as suggested 
by the observation of glitches in the rotation period [73]. We do not have a model for 
this variation and, therefore, it is not currently possible to evaluate the photon spectra 
at the present day, but we can use the Crab photon spectra observed after ~ 1000 years 
from the birth. For the Crab-like source the initial spin-down time scale is tq ~ 700 
years and the characteristic age is Tage,c ~ 11 kyrs. The conversion efficiency of the 
spin down power into particle luminosity ( 77 ) is defined by {ll}{rie*+rie- - = 0 . 1 ). 

Therefore, we take electron-positron spectra, normalized to the photon emission like in 
the Crab nebula, and propagate the source spectrum after T ~ 1000 years, as assumed 
in [18]. In Fig. 3 and 4 the positron or electron spectrum, obtained from Eq. (13) 
using all the Crab parameters, is shown for Vela-X (Model 2). The main parameters of 
the two models are summarized in Tab. 3. 

From an inspection of Fig. 3, one may remark that the measured electron and 
positron intensities can be accounted for by the flux expected from Vela-X within mod¬ 
els 1 and 2. Fig. 4 shows the comparison between our models (cLIS and Vela-X con¬ 
tributions) and the absolute AMS-02 electron plus positron flux. The shadowed bands 
keep into account the uncertainties due to the cLIS parameters (see Tab. 2), the error 
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Figure 3: Two model of this analysis of the expected positron or electron omnidirectional intensities from 
Vela-X compared with the half electron plus positron signal as in Fig. 1. The shadowed bands reflect the 
uncertainty on Vela distance, that is about 6%, and the variation due to the assumptions on the initial rotation 
period. 

on Vela-X distance and, for Model 1, also the variation due to the assumptions on the 
initial rotation period. In the top graph, the band is in agreement with the AMS-02 data, 
thus, Vela-X alone can account for the electron and positron excess components. In the 
bottom one, for energies between 20 and 80 GeV, the AMS-02 data are slightly above 
the model. This discrepancy can be accounted for introducing aged pulsars contribut¬ 
ing with lower energy electrons and positrons (in comparison with Vela-X), while the 
photon emission of their PWN is no more observable. The main contribution coming 
from these mature pulsars is due to Monogem. Its electron and positron contributions, 
evaluated with Model 2 (Monogem treated as Crab-like), is a factor 2 higher than the 
Vela-X one at 40 GeV. 
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Figure 4: Electron plus positron contribution from cLIS and Vela-X Model 1 and 2 (solid lines) compared 
with the AMS-02 flux. The bands keep into account the uncertainties of the cLIS parameters (see Tab. 2), 
the error on the Vela-X distance and, for Model 1, also the variation due to the assumptions on the initial 
rotation period. 




5. Discussion and Conclusions 


In the present paper, we analysed the AMS-02 electron, positron and electron plus 
positron spectra in the energy range above ~ 10 GeV. The excess spectra for electrons 
and positrons were obtained by subtracting the expected “classical” LIS’s, computed 
with GALPROP, to the omnidirectional distributions observed by AMS-02. A compar¬ 
ison above ~ 50 GeV indicates for electrons and positrons the same flux and the same 
slope. These excess spectra can be accounted for by pulsar sources in which generated 
electron-positron pairs will be accelerated by the surrounding pulsar wind nebula. In 
particular, we evaluated electron and positron spectra generated in the Vela-X PWN 
and propagated them to the Earth. We used two different models built using observed 
parameters of Vela and Crab nebulae. Finally, we compared results with observations. 
Both models, taking into account uncertainties and assumptions, are not in disagree¬ 
ment with the AMS-02 excess components at energy larger than ~ 100 GeV. Vela-X 
Model 2, built on Crab, requires a particle conversion efficiency which is an order 
of magnitude higher than Model 1. Our results are also in agreement with a smooth 
change in the spectral shape as reported by AMS-02. 

Comparison with data, for energy above 100 GeV (see solid lines in Fig. 4), indi¬ 
cates the possibility of an extra source similar to Vela. There is only another known 
pulsar, B1737-30 (or J1740-3015), with parameters similar to Vela. It is 400 pc far and 
20600 years old [75, 76], but the photon emission of its nebula has not been observed 
yet. Moreover, statistically we expect just one (or very few) more pulsar like Vela: we 
can consider a pulsar birth rate of 0.9-2 objects per century [77, 78], a spatial distribu¬ 
tion like in [77] and [79], pulsars with age from 10 to 50 kyrs in a volume of about 1 
kpc^ around the Earth. 

At energies lower than 100 GeV, the tiny discrepancy of Fig. 4 (Model 2) can be 
covered by mature pulsars (10^ years) because the low energy electrons and positrons 
released are still diffusing to the Solar System. The main contribution from this sources 
comes from Monogem and it is a factor 2 higher than the Vela-X one at 40 GeV. 

Our models predict that a single close PWN, Vela-X, may be responsible for at 
least half of the electron and positron excess in CRs, thus, a dipole signal in the CR 


17 



arrival directions could be detected. If no other source is missed, we expect a dipole 
anisotropy above 100 GeV centered in the direction of Vela. At 200 GeV, the dipole 
anisotropy for electrons plus positrons from Vela-X is expected to be of the order of 
~ 2% [80] not yet excluded by the FERMI experiment [81]. Conversely, at lower en¬ 
ergy, several sources can contribute to the electron and positron spectra, but the angular 
distribution of all these sources should be more isotropic. AMS-02 observation could 
be able to provide additional data to the anisotropy studies. 

Our models predict a change of the slope at energies below 10 TeV (see Fig. 3). 
Future measurements up to this energy range can be used to confirm or reject such a 
prediction. 

Our results are in agreement with models describing the origin of the pulsar wind 
nebula (e.g., the one discussed in 40), where the minimal electron and positron lumi¬ 
nosity is not expected to be lower than a few percent. Thus, in this context, the current 
results can be used to constrain the fraction of the spin-down luminosity which is trans¬ 
ferred to particle acceleration needed to fit the excess spectra observed by AMS-02. 

Finally, the satisfactory agreement between the models and the data leads to keep 
into account PWN as source of electrons and positrons. Therefore, a realistic FIS 
should include this type of electrons and positrons sources. 
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